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Abstract 

The stability of the 1+1 dimensional solution of Israel- Stewart theory is in- 
vestigated. Firstly, the evolution of the temperature and the ratio of the 
bulk pressure over the equilibrium pressure of the background is explored. 
Then the stability with linear perturbations is studied by using the Lyapunov 
direct method. It shows that the shear viscosity may weaken the instabil- 
ity induced by the large peak of bulk viscosity around the phase transition 
temperature T^. 
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1. Introduction 



A hot and dense partonic matter has been created at the Relativistic 
Heavy Ion Colhder (RHIC) of Brookhaven National Laboratory which is 
thought to be a kind of perfect hquid of quark and gluon i2|. The study on 
the properties of this matter is a hot issue. Presently one property which is of 
very interest is the exact value of the ratio of viscosity over entropy density 
for the matter. To extract the value, one needs to cornpare experimental 
data with relativistic viscous hydrodynamics simulation |3|, 14, [sl, la l7|. 

Two different kinds of relativistic viscous hydrodynamics have been devel- 
oped so far. One is the so-called first order relativistic viscous hydrodynamics 
which was first developed by Eckart and varied by Landau and Lifshitzjs, 9|, 



10|. But there are two 



and many calculations have been done since then 
problems in the approaches: one is that dissipative fluctuation may propa- 
gate at a speed larger than the speed of light and thus leads to a causality 



problem; the other is that the solution which may develop instabilities 
These problems have been studied extensively in Refs. [3| • The other 
is the second order theories among which the Israel-Stewart theory is pop- 
ularly used so far jl^. However, before Israel-Stewart theory is applied to 
describe heavy ion collisions, one should know whether it is stable or not 
Till now, there have been some stability analysis on the issue 
In jisl, Q, 0, 18 1, the authors used the plane wave perturbation method to 
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16 



17 



18|- 



study the stability of the theory around a hydrostatic state and discussed the 
regime of validity of hydrodynamics. Also in using the Lyapunov direct 



method as in 



13|, the authors analyzed the stability of the scaling solution 



with only bulk viscosity and presented the stable regions of the scaling solu- 



2 



tion for both homogeneous and inhomogeneous perturbations. And several 
groups have compared this theory with the experimental data and show that 
this theory can be used without aproblem when the shear viscosity and bulk 



viscosity are small 



SHE 



161,13, 



19|. 



It is well known that hydrodynamics fails to describe the HBT results, 
i.e, the calculated ratio of HBT radii in outward direction over that in side- 
ward direction is higher than the experimental data 20| . Recent att ernp t to 
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reconcile them is to take bulk viscosity into account. But in 
the authors found that the HBT radii hardly change if the peak of bulk 
viscosity is small and the other parameters such as the equation of state, 
initial condition and so on keep unchanged. However, in 23| the authors 
analyzed the stability of Navier-Stokes theory with the bulk viscosity which 
has a large peak around Tc and found that there exists some inhomogeneous 
modes which will tear the system into droplets. If this is the case, it will help 
us resolving the interferometric data 2J]. But there are some basic problems 
with Navier-Stokes theory as mentioned above and the authors did not con- 
sider the effects of shear viscosity additionally. In this paper we use the same 
method to study if this could happen in the Israel- Stewart theory and what 
is the role of the shear viscosity. 



2. Israel-Stewart theory and linear perturbations 

The general hydrodynamic equations arise from the local conservation of 
energy and momentum j25| 



0, 



(1) 



3 



where the energy-momentum tensor without heat conduction is decomposed 



into the following form 

T^'^ = eu^M'^ - (p+n)A'^^ + vr'^^ (2) 

where e and p are the local energy density and thermal equilibrium pressure, 
and is the 4- velocity of the energy flow which obeys u^u^ = 1. 11 is the 
bulk viscous pressure, and A^'^ = g^^'^ — u^u'^ is transverse to the flow velocity, 
that is A'^'^Uu = 0. vr^'^ is the traceless shear viscous pressure tensor. With 
Eq. ([1]) and Eq. ([2]) , we can get the evolution equations of the energy density 
and the 4-velocity of energy flow. 

In Israel- Stewart approach, the kinetic evolution equations of the bulk 
pressure 11 and the traceless shear viscous tensor n^" are j3] 

DU = -l(n + CV-n), (3) 



where D is the time derivative in the local fluid rest frame and fulfills D = 



df^. The angular bracket notation is defined by {V^u^) = ^{V^u'^+V^u^) — 
|(V-n)A^'^. 1] and ( denote the bulk and shear viscous coefficient, respec- 
tively, and Tn are the relaxation time for the bulk pressure and the shear 
tensor, respectively. They can be related to r] and ( as follows [ll| 



= 277/32 , Tn = C/^o- (5) 

where P2 and /3o are the relaxation coefficients that need to be calculated 
from other theories. 
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When hydrodynamics is apphed to describe heavy ion colhsions, people 
always use the symmetry of the collision system to simplify the evolution 



equations. T 
2+1 solution 



lus we have the scaling solution, the 1+1 solution and the 
Here we focus on the 1+1 solution of the Israel- 



Stewart theory. For 1+1 solution, the system has boost-invariant longitudinal 
expansion and transverse expansion in one dimension. For these geometries, 
it is convenient to work in the co-moving and radial coordinates r, r, 0, r]. 
We consider a small perturbation, and the backgrounds eo(T), 7r(^Q(r), no(T) 
evolve with time 

e(r,r) = eo(r) + (5e(r,r), (6) 
v{t, r) = 6v{t, r), (7) 
KiT,r) = 7r%iT)+5n^^iT,r), (8) 

U{T,r)=Uo{T) + 6U{T,r). (9) 

With Eqs. ([IHl]) and Eqs. ([HHHD the evolution equations of backgrounds 
can be given as follows 

1 1 

9reo = -(eo+Po + no)- + 7r:jo-, (10) 

r r 



T^o^r^o + <o = (11) 



TnodrKo + Ko = 'l — , (12) 

o r 
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rnoa.no + Ho = (13) 
r 

We can also get the evolution equations of the perturbations 

[(eo + po + no - -nlojdr + dr{pQ + Ho) - {dr + -Xo + -Ko\^^ 

T T ' 

+ dri6p + 511) - {dr + ^)5< - = 0, (14) 

[(eo +P0 + Ilo){dr + -)- Koidr ' -) + -<o]^^ + drSe 

+ -(5e + (5p + (5n)-i(5< = 0, (15) 
r r ' 

r.o9.K + 5< = _^(4!^ _ ) _ 2 „(a.5^; + ^) + 1^, (16) 
' ' r^ro o r d r 6 T 

rMi^r + (^TT^ = (g — + ^ro) - i:VQ[-^Or5v + — ) - - — , (17) 

r,rO <J T <J r o r 

rno^.OT + = ^(Ho + ^) - Co(5.5t; + -) - ^. (18) 
Tno T r r 

In order to get rid of the space-like derivatives, we can do the Hankel trans- 
form (Fourier-Bessel transform) due to the geometry system we use hereQ]. 
After introducing 

Sf, = {dr + ^)5< + U'k'^ (19) 

we do the following Hankel transforms, 

5v{T,r)= dkJi{kr)k6v{T,k), (20) 
Jo 



oo 



Se{T,r)= dkJok{kr)6i{T,k), (21) 







/■oo 

Sn{T,r)= dkJi{kr)k6n{T,k), (22) 
Jo 



STrV{T,r)= I dkMkr)k6n^{T,k), 
'o 



(23) 



5n(r,r) 



dkjQ{kr)k5Il{T, k), 



(24) 



where k is the wave number. It stands for a homogeneous perturbation when 
k equals zero or an inhomogeneous perturbation when k is not zero. Its 
range for a realistic QGP fluid has been roughly estimated in 
Eqs. dHHM]), we can get the evolution equations for 6v, 6e, 611^, 671, 611 



13|. With 



[(!eo +po + no - <o)'9r + dr{po + no) - {dr + -)Ko + -Ko\^^ 

T T 



kclde - kdli - 5^ = 0, (25) 



(eo + Po + Ho + ^^o)^^ + {dr + ^^)5e + -fi - U^"^^ = 0, (26) 



2^ 

3 TVo 



1 Am 



4 1 d7]^ 
3rr^o 5e 



+ {dr + —)6n^, = 0, (27) 



3 r^o 



, r~ 2 1 .ar], , 

iofcae+ [ — )Qkde 

3rr^o <9e 

+ {dr + — )5^ = 0, (28) 



- ^(^ + no)(^)o<5e + ^(^)o^e + (9. + — )5II = 0. (29) 
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3. Results 



In order to get numerical results, the values of viscosities and the relax- 
ation time are set up as described below. As to shear viscosity, the strong 
coupling theory and the hydro dynamic and transport model show that - can 



not be too large 



27 



28|. Here we use two different values ^ = 0.02 and 



0.2 to see its effects. The relaxation time of the shear viscosity is set to the 
Boltzmann gas result = in our work. Although recent results of SU(3) 
Yang-Mills theory show that the effects of a potentially large bulk viscosity 
near Tc are more subtle to detect in spectral integrals 29| than previous work 
suggested 301, but there are still some possibilities that - becomes large 

n n ' 

|31| . So the following parametrization as in [23| is adopted to see 



around Tc 

the effects of bulk viscosity on the stability 



pQCD 



+ 



=exp[ 



2a2 



(30) 



V2na 

where t = T — Tc, cr = O.OlTc, ZpQcn ~ 10^^ and different zq denotes 
different magnitudes of (/s. For the relaxation time of bulk viscosity, the 
parametrization similar to the relaxation time of shear viscosity is employed 

c 



Tn 



e + p 



(31) 



where different b is used to see the effect of relaxation time of bulk viscosity. 

To connect to relativistic heavy-ion collisions, the equation of state in 32 1 
is adopted and the initial temperature Tq is set to be 0.34 GeV. Firstly, let 
us see the effect of viscosity on the evolution of background. Fig 1. shows 
the profiles of background temperature for different rj/s and for different zq 
and b, corresponding to different magnitudes and relaxation times of bulk 
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viscosity. Each plot consists of two different rj/s. From the left to the right, 
the plots represent zq = 0.0 and 0.1 Tc, respectively. These plots show that 
the shear viscosity can slow down the evolution of the source. We can also see 
that when b is large which means that the relaxation time of bulk viscosity 
is large, bulk viscosity hardly have effects on the evolution of temperature 
even when | has a large peak. But when the relaxation time is small, for 
a larger magnitude of bulk viscosity, we can see that the system stays at a 
nearly constant temperature which is the same as the case found in Navier- 



Stokes theory in 



231]. This effect is not hard to be understood. Because the 



role of the relaxation time is to delay the appearance of viscous forces 
bulk viscosity has effects on the evolution of the source only after a time- 
scale Tu- The larger the relaxation time is, the longer time the bulk viscosity 
needs to affect on the evolution of the system. Meanwhile, the system is still 
evolving and it may be cooled down to be below Tc where bulk viscosity can 
be negligible. So if the relaxation time is large enough and the bulk viscosity 
has a small width, bulk viscosity may have negligible effects on the evolution 



of background. It is consistent with the results in 22] . In [22!] we found that 
the width of bulk viscosity has larger effects on the evolution of the source 
than the magnitude of bulk viscosity. 

Fig 2. shows the profiles of the ratio of bulk pressure over the equilibrium 
pressure of the background for different cases as in Fig 1. We can see that 
U/p has a peak when (/s has a peak around Tf.. And the fact that U/p 
exceeds one means that the state is far away from equilibrium. It indicates 
that the matter is not only hydrodynamically unstable, but also thermody- 
namically unstable as stated in when the peak of bulk viscosity is large 
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and the relaxation time of bulk viscosity is small. We can also see that shear 
viscosity will decrease II/p which indicates that shear viscosity will weaken 
the instability that induced by bulk viscosity. This behavior is also found in 
the next analysis about linear perturbations. 

Now we use Lyapunov direct method to study the stability of 1+1 so- 
lution of Israel- Stewart theory. The evolution equations of perturbations 
Eqs. (125K I29|) can be rewritten in the following matrix formula 

dJY = A6Y, (32) 

where 

^ 6v ^ 

6i 

SY= Snll . (33) 

When one uses the Lyapunov direct method, the Lyapunov function should 
be given first. The Lyapunov function must be positive definite. If it is 
a monotonically decreasing function, then the solution is stable; if it is a 
monotonically increasing function, then the solution is unstable. The more 

n □ 

detailed description about this method can be found in[13|, [ITJ. Here we 
assume the Lyapunov function is = 6Y'^6Y. Then the evolution equation 
of the Lyapunov function V is 

d,V = SY'^iA^ + A)5Y, (34) 



after a short derivation, we can get 



(35) 



10 



where Amax and A^m are the largest and smallest eigenvalues of A + , 
respectively. It can be shown that the solution is stable when Amax < and 
is unstable when Amin > 0. Fig. 3 and Fig. 4 show the values of Amax and 
Amin for k = and k = 3, respectively. Each figure has differen cases as 
stated for Fig. 1. We can see that neither the stable region nor the unstable 
region can be determined because Amax is always larger than zero and Amin 
is always smaller than zero. As for the scaling solution, the unstable regions 
also can not be found but the stable regions can be determined |l7|. These 
plots show that Amax and Amin have rapid change when the bulk viscosity is 
large and the relaxation time is small. We can also see that shear viscosity 
will delay this kind of rapid change and reduce the magnitude of Amax and 
Amin- This phenomenon is the same for different k, but the magnitude of 
Amax and Amin is larger with larger k. 

We can see from Fig. 3 and Fig. 4 that the absolute value of Amax 
and Amin increase rapidly which means the growing and damping rates of 
perturbations increase rapidly due to the large peak and small relaxation 
time of bulk viscosity. Therefore the perturbations may rapidly grow to a 
value comparable with the background. So they will break local homogeneity 
and play an important role in the subsequent evolution of the system. The 
created inhomogeneities have no reinteractions. It is possible that isolated 
fragments will be created and move away from each other. In 23(|, the authors 
argued that this may be a reason that the source will be clusterized and then 
decoupled as the fireballs. By adding a further free parameter which is the 



cluster size to the system, this mechanism may solve the HBT problem 



24|. 



In the Israel-Stewart theory this phenomenon may also happen when the 
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peak of C/s is large, rn is small and C/'^ is not too large. 
4. Conclusions 

In summary, the stability problem of the 1+1 solution of Israel- Stewart 
theory is studied. Firstly, the evolution of the temperature and the ratio of 
the bulk pressure over the equihbrium pressure are studied. We find that 
both shear and bulk viscosity slow down the evolution of temperature. And 
the ratio of bulk pressure over the equilibrium pressure (II/p) will exceed one 
with a large peak of bulk viscosity, which indicates the state is unstable. The 
shear viscosity reduces the magnitude of II/p to weaken the effects of bulk 
viscosity. Then using Lyapunov direct method, we can not determine the 
stable or unstable regions. We also find the phenomenon which may drive 
the source to clusterize which is similar to that in Navier-Stokes theory. 
However, this phenomenon will happen only when the peak of bulk viscosity 
is large enough, the relaxation time of the bulk viscosity is small and the 
shear viscosity is not too large. 
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Figure 1: Evolutions of background temperature for different shear and 
bulk viscosity sets. See texts for details. 

Figure 2: 11 /p as a function of T/Tg for different shear and bulk viscosity 
sets. See texts for details. 

Figure 3: Xmax and Xmin ior k — with different shear and bulk viscosity 
sets. See texts for details. 

Figure 4: Xmax and Xmin ior k — 3 with different shear and bulk viscosity 
sets. See texts for details. 
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Figure 1: Evolutions of background temperature for different shear and bulk viscosity sets. 
See texts for details. 




Figure 2: li/p as a function of T/Tc for different shear and bulk viscosity sets. See texts 
for details. 
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Figure 3: Amax and Xaun for k = with different shear and bulk viscosity sets. See texts 
for details. 



4000 - 



2000 



-2000 



' I ' I ' I 




S = =-=^ ... ^ Bi W W 



6000 



I ' I 

Zo=0.1T 



0.8 0.9 1 1.1 
T/T 




ri/s=0.02 b=0.2 
ri/s=0.02 b=6.0 
ri/s=0.2 b=0.2 
ri/s=0.2 b=6.0 



Figure 4: Amax and Amin for A; = 3 with different shear and bulk viscosity sets. See texts 
for details. 
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